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Abstract — The paper studies the inferences of wave equations for electromagnetic fields when 
there are gravitational fields at the same time. In the description with the algebra of octonions, 
the inferences of wave equations are identical with that in conventional electromagnetic theory 
with vector terminology. By means of the octonion exponential function, we can draw out that 
the electromagnetic waves are transverse waves in a vacuum, and rephrase the law of reflection, 
Snell's law, Fresnel formula, and total internal reflection etc. The study claims that the theoretical 
results of wave equations for electromagnetic strength keep unchanged in the case for coexistence 
of gravitational and electromagnetic fields. Meanwhile the electric and magnetic components of 
electromagnetic waves can not be determined simultaneously in electromagnetic fields. 



1. INTRODUCTION 

In the electromagnetic field, the inferences of wave equation are believed to be correct. Until to now, 
this validity is confined only to one special condition that there exist only the electromagnetic fields. 
When there are gravitational fields, the people doubt whether the wave equation of electromagnetic 
field is still correct or not. The validity about wave equations is being questioned continuously, and 
it remains as puzzling as ever. The existing theories do not explain why the deductions of wave 
equation should keep unchanged, and then do not offer compelling reason for the unique situation. 
The paper attempts to find out why the deductions of wave equation keep the same in most cases, 
even in the gravitational field. 

The law of reflection (or specular reflection) is the mirror-like reflection of light from a surface, 
and the angle of incidence equals the angle of reflection [IJ [2] . This behavior was first discovered 
through careful observation and measurement by Heron of Alexandria, Greek. The law of reflection 
was dealt with the Principle of the Shortest Path of Light, which was later stated in this form by 
P. Fermat in 1662 [3]. The most modern form is that the path is at an extremum [H [5]. 

Snell's law (also known as law of refraction) was first accurately described in a mathematical 
form by Ibn Sahl of Baghdad. It was rediscovered by T. Harriot in 1602, who did not publish his 
results. In 1621, W. Snell derived a mathematically equivalent form, that remained unpublished 
during his lifetime, until C. Huygens mentioned Snell's discovery in 1678 0E]- 

A.- J. Fresnel rediscovered interference and provided a mathematical theory of diffraction based 
on the wave theory [HJ [9]. In 1817, T. Young proposed a small transverse component to light, while 
yet retaining a far larger longitudinal component. By the year 1821, Fresnel was able to show via 
mathematical methods that polarization could be explained only if light was entirely transverse, 
with no longitudinal vibration whatsoever. He relied on mathematical analysis which enabled him 
to make predictions which could be accurately verified and he removed many of the objections to 
the wave theory of light [TUl E] • 

The algebra of quaternions was invented by W. R. Hamilton in 1843 |12j . and then was first used 
by J. C. Maxwell to represent field equations of electromagnetic field in 1861 [13]. O. Heaviside 
in 1884 recast Maxwell equation in terms of vector terminology and electromagnetic forces |14j . 
thereby reduced the original twenty equations down to the four differential equations. 

In 1871, H. Helmholtz made clear the electromagnetic theory effectively [15(. I16j . And Helmholtz 
equation and boundary conditions can deduce some wave features, including the law of reflection, 
Snell's law, Fresnel formula, etc |17j . Later H. R. Hertz became famous as the first to demonstrate 
the electromagnetic radiation. 

In the paper, in terms of the quaternion feature, we obtain field equation and related wave 
feature in the gravitational field. Rephrasing with the algebra of octonions, we achieve Maxwell 
equation in electromagnetic fields, and then some wave properties including the wave equation, law 
of reflection, Snell's law, Fresnel formula, and total internal reflection, etc. 
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2. WAVE EQUATION FOR FIELD STRENGTH IN GRAVITATIONAL FIELD 

The wave feature of gravitational field can be described by the algebra of quaternions. In the 
quaternion space, the coordinates are ro, r\, and rs, with the basis vector K g = (1, i\, i%, 83). 
The radius vector is M g = S(rjtj), and the velocity Y g = S(^jij), with = 1. Herein ro = v$t . v$ 
is the speed of gravitational intermediate boson, and t is the time, i = 0, 1, 2, 3. j = 1, 2, 3. 

In the quaternion space for the gravitational field, the gravitational potential is, A g = S(aj«j). 
And the gravitational strength is M g = = (} o A g .The gravitational strength M g covers two 

components, g/^o = <9oa + Va > ana ^ b = V x a . The gauge equation is &o = <9oao + = ^> with 
a = Yi(djij) and V = Y>(ijdj) . Herein the o denotes the quaternion multiplication, the operator 
= T,(iidi), with di = d/dri . 

Comparing with classical gravitational theory, the linear momentum density, P = mV g , is the 
source S g of gravitational field. And the latter one is defined as follows, 

(Mg/v Q + O M g = -flgSg , (1) 

where m is the mass density. * denotes the quaternion conjugate. [x g = A-kG/vq is a coefficient, 
and G is the gravitational constant. IB* oM g / (2/i ff ) is the gravitational field energy density, with 
the extra mass part Am = B* o M g /(vQ/j, g ). 
Further, the above can be separated into 

V • b = , d b + V* x g/v = , (2) 
V* • g = ~mj e g , <9og/wo + V* x b = -/i g s , (3) 

where the coefficient e g = l/(/j, g VQ). s = T,(sjij). Sj = pj = rnvj. so = Po = fhvQ. 711 = 171 + Am. 

The above means that the Newton's law of gravitation in the quaternion space |18| is the same 
as that in classical gravitational theory \19\ [20] . 



Table 1: The quaternion multiplication table. 
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2.1. Quaternion exponential function 

In the quaternion space, defining the complex quaternion as 

exp {i(a + ib)} = cos{a + ib) + isin(a + ib) , (4) 

where a and b are all real, i is the imaginary unit, with i 2 = — 1 . Meanwhile the quaternion unit 
i is one function of ij , with i o i = — 1 . 

When 6 = 0, the above reduces to the quaternion, exp(ia) = cos(a) + isin(a) . When a = , 
Eq.(4) reduces to the Hamiltonian bi-quaternion, 

exp(iib) = cos(ib) + \sin(ib) . (5) 

The bi-quaternion exp(-iia) satisfies, 

exp(-iia) o exp{i\a) = 1 , exp(-iia) o exp(-iia) = exp(—ii2a) , (6) 
djexp(-iia) = iikj o exp(-iia) , d 2 exp{—i\a) = kjexp(-iia) , (7) 

where a = —T,(kjrj), with kj being the real coefficient. 

In case of neglecting the direction of i , Eq.(4) reduces to the complex exponential function, 

exp{i(a + ib)} = cos(a + ib) + isin(a + ib) = exp(ia)exp(—b) . (8) 
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Table 2: The operator and multiplication of the physical quantity in the quaternion space. 



definition 


meaning 


Va 


-(diai + d 2 a 2 + d 3 a 3 ) 


V x a 


ii(d 2 a 3 - d 3 a 2 ) + ia{d 3 ai - #ia 3 ) + i 3 (<9ia 2 - <9 2 ai) 
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iid ai + i 2 d a 2 + i 3 d a 3 



2.2. Wave equation 

Researching the transmission problem of gravitational waves needs to solve field equation. In 
a vacuum far away from the gravitational sources, there does not exist the mass so and linear 
momentum s, and then field equation can be reduced to, 

0* o M g = . (9) 

Applying the operator to the above, we get the wave equation 

o (0* o Mg) = , or {dl + d\ + d\ + dl)M g = . (10) 

The above is the Laplace equation in the Euclid space, meanwhile it is the wave equation in 
the quaternion space. And this Laplace equation can also be obtain from Eq.(9) by means of the 
conjugate operator 0* and the field equation. 

Proceeding with the operator do and V, we can obtain the wave equation about the components 
of field strength from the field equation Eqs.(2) and (3) directly, 

(d 2 + d 2 + d 2 + d 2 ) g = o, (d 2 + d 2 + dl + d!)h = o. (ii) 

The above means that Eq.(10) can be separated into two parts in Eq.(ll) when §> g = 0. And two 
strength components of gravitational field, b and g, are both possessed of wave features, although 
their wave impact may be too tiny to be detected at present. 

2.3. Transverse wave 

For the gravitational wave with the angular frequency u>, the field strength should be a harmonic 
function, coscot. In order to operate handily and cover more physics contents, it is convenient to 
substitute the function exp(-iojt) for the cosojt. 

The gravitational strength g and b can be written as follows, 

g = g(r)exp(—iujt) , b = h(r)exp(— iujt) , (12) 

combining this with Eq.(ll), gives 

{-(u/vo) 2 + Xdf } g(r) = , {-(uj/v ) 2 + Sdf } b(r) = . (13) 

From the analysis of field equation, we find that the amplitude of gravitational strength will be 
increased steadily or deceased continuously. And then the field strength should be one hyperbolic 
cosine, cos(ia). It will substitute the function exp(-iia) for the cos(ia). Herein the wave vector 
k = Y,(ijkj), while the vector radius r = T,(ijrj), and ia is an imaginary angle. 

And then, the gravitational strength g(r) and b(r) are 

g(r) = g o exp(-iia) , b(r) = bo o exp(-iia) , (14) 

where g and b both are constant vectors in the quaternion space. 
Further, substituting the above in Eq.(13), we have the result, 

- {u/vo) 2 + m 2 = . (15) 

In the gravitational field, we will find that gravitational waves are the transverse waves in a 
vacuum. From the field equations 



V-g = , V-b = , 



(16) 
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we obtain 

k • g = , k • g' = , k • b = , k • b' = , (17) 

where g' = g o i , b = b o i . 

The above states that the gravitational wave components, g and bo , belong to the transverse 
waves. So do these two new wave components, g and b . Moreover, the amplitudes of wave 
components, g and b' , are the same as that of g and bo respectively. 

In the gravitational field, there are two new kinds of gravitational waves. The field equations, 

V xg = dh/dt , V-g = , (18) 

can be rewritten as, 

Vog = db/dt. (19) 

Expanding the exp(-iia), the above implies that 

[k o g + ujho]cosh(a) + [k o g — ujh' ]sinh(a) = , (20) 

considering Eq.(17) and the a can be an arbitrary value, and then 

k x g + ub Q = , k x g - wb = , (21) 

where sin(ia) = isinh(a) , cos(ia) = cosh(a) . 
In the same way, from the field equations 

ugV x b = dg/dt , V • b = , (22) 

we get 

k x b' + Ld(g /v )/v = , k x b - uj(g' /v )/v = . (23) 

The above means that there are some relationships between M g with M g o i . The k and g will 
yield a new component b , while the k and bo produce the g' . The g and b' can be determined 
at the same time, while the bo and g can be measured simultaneously. Provided there are some 
substances in the nature may slow down the speed of gravitational waves, the field strength will 
possess the wave features, including the reflection, refraction, dispersion and energy flux etc. 



Table 3: The octonion multiplication table. 
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3. WAVE EQUATION FOR FIELD STRENGTH IN ELECTROMAGNETIC FIELD 

The wave feature of gravitational field and electromagnetic field can be described simultaneously 
by the octonion space, which consists of two quaternion spaces. 

In the quaternion space for the electromagnetic field, the basis vector is E e = (Jo, Jj, I2, I3), 
the radius vector is R e = S(i?jJj), and the velocity is V e = S(ViJ). The E e is independent of the 
E 9 , with E e = E g o I and Jq o J = — 1. Vq is the speed of electromagnetic intermediate boson. 

These two quaternion spaces can be combined together to become an octonion space [5TJ . with 
the octonion basis vector E = (1, «i, 13, Jq, Ii, I2, 13). The radius vector in the octonion space 
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is R = S(rjij + Rili), and that the octonion velocity is V = + V^Jj) . Herein the symbol o 

denotes the octonion multiplication. 

The potential of gravitational field and electromagnetic field are A g = £(ajij) and A e = £(Aj Jj) 
respectively. They are combined together to become the field potential, A = A g + k eg A e , in the 
octonion space. The field strength IB consists of the gravitational strength M g and electromagnetic 
strength B e . That is,B = <0°A = B 9 + k eg M e , with the k eg being a coefficient. The strength B e 
covers two parts, E/Vb = (doA + VAo) ° Io , and B = (V x A) o I . The gauge equations are 
b = and B = 0. Herein B = d A + V • A, with A = 

Comparing with the classical electromagnetic theory, the electromagnetic source is the current 
density § e = q¥ e , while the gravitational source is the linear momentum density. And the source 
§ was devised to describe consistently the field source of electromagnetism and gravitation. 

(M/v + 0)* o B = -fjS = -(jigSg + k eg n e S e ) , (24) 

where fi is a coefficient, and fj, e is the electromagnetic constant, with k 2 eg = n g /fJ- e ■ Q is the electric 
charge, and * denotes the conjugate of octonion. 

According to the basis vectors, the above can be decomposed further as follows, 

B* oM/v + 0* °B 9 = -Hg^g , 0* o B e = - (l e S e , (25) 

where the former equation is similar to Eq.(l) and is suitable for the gravitational field, while the 
latter equation is for the electromagnetic field. 

Further, the latter equation in the above can be rewritten as follows 

V • B = , a B + V* x E/Vfc = , (26) 
V* • E = -(q/e e )I , d B/V + V* x B = -^S , (27) 

where the coefficient e e = l/( / u e V 2 ). S = T>(SjIj), Sj = qVj. So = Solo, So = qVo- 

By contrast with Maxwell equation in terms of vector terminology, it is found that the above is 
the same as that in the conventional electromagnetic theory respectively, except for the direction 
of displacement current. And the gauge equation Bq = is different. 



Table 4: The quaternion, Squatcrnion, and octonion with the complex number. 

subspace meaning coordinate basis vector 

quaternion quaternion real number i$ , i\ , ia , i% 

bi-quatcrnion imaginary number «o, «i,i2,*3 

complex quaternion complex number «o j H > h , *3 

Squaternion Squaternion real number Jo, Ji, l2,Is 

bi-Squaternion imaginary number Jo, I\, 12,13 

complex Squaternion complex number Jo, Jl, J2, J3 

octonion octonion real number io , ii , ia , is , Jo , Jl , I2 , 13 

bi-octonion imaginary number io , i\ , h , h , Jo , Ji , J2 , J3 

complex octonion complex number io , i\ , h , is , Jo , Jl , h , I3 



3.1. Octonion exponential function 

Similarly the octonion O is called as the bi-octonion when the coordinate is the imaginary number. 
Meanwhile the octonion O is the complex octonion with the complex number coordinate. 
The octonion O is separated into two parts, the quaternion Oi and Squaternion O2 , 

O = <0>i + 2 , 2 = Io o O3 , (28) 

where O3 is the quaternion. Squaternion means the second quaternion or step-quaternion. 
In the octonion space, defining the complex Squaternion as 

exp{I(a + ib)} = Io o {cos(a + ib) + i'sin(a + = Iocos(a + ib) + Isin(a + ib) , (29) 

where a and b are all real. The i is an imaginary unit, with i 2 = — 1 . The quaternion unit i' is one 
function of ij , with i' o i' = — 1 . Meanwhile the unit I = Iq o i' , with I o I = — I . 
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When 6 = 0, the above reduces to the Squaternion, exp(Ia) = Iocos(a) + Isin(a) . Similar to 
Hamiltonian bi-quaternion, when a = 0, Eq.(29) reduces to the bi-Squaternion, 

exp(Iib) = Iocos(ib) + Isin(ib) , (30) 

with 

exp(-Iib) = Iocos(ib) — lsin(ib) . (31) 

The bi-Squaternion exp(-ila) satisfies, 

exp{— ila) o exp(ila) = io ° exp(iI2a) , exp(ila) o exp(ila) = — 1 , (32) 
djexp(-ila) = HKj o exp(-ila) , d 2 exp(—ila) = K 2 exp{—i\a) , (33) 

where a = —T,(Kjrj), with Kj being the coefficient. 

In case of neglecting the direction of io ; the imaginary unit i will replace the unit Iq , and then 
Eq.(29) reduces to the quaternion exponential function similar to Eq.(4). 

Table 5: The operator and multiplication of the physical quantity in the octonion space. 



definition 


meaning 


VS 


-(diSi + d 2 s 2 + d3S 3 )Jo 


V x S 


-Ii(ftS 3 - d 3 S 2 ) - 7 2 (a 3 5! - d lt S 3 ) - - &Si) 


VSo 


Ji9i5 + hd 2 S + I 3 d 3 S 


d S 


hdoS! + hd s 2 + i 3 d s 3 



3.2. Wave equation 

Studying the wave transmission feature in electromagnetic fields needs to solve Maxwell equation. 
In a vacuum far away from the electromagnetic sources, there does not exist free electric charge 
density So and current density S, and then field equation can be reduced to, 

0* o B e = . (34) 

Applying the operator to the above, we have Laplace equation, 

o (0* o B e ) = , or (dl + d\ + dl + <9f)B e = . (35) 

The Laplace equation is the wave equation in the octonion space. This wave equation can be 
obtain from Eq.(34) with the conjugate operator ()* and the field equation also. 

Proceeding with the operator do and V, we can obtain the wave equation about the components 
of field strength from Maxwell equation Eqs.(26) and (27) directly, 

(d 2 + d 2 1+ d 2 2 + d!)B = o , (d 2 + d 2 + d 2 + d 2 )B = 0. (36) 

The above means that the Eq.(35) can be decomposed into two equations in the above. And 
two strength components of electromagnetic field, E and B, are both possessed of wave features, 
and can be detected at present. 

3.3. Transverse wave 

In a similar way, for the electromagnetic wave with the angular frequency u, the field strength 
should be a harmonic function, cosoot, and can be chosen as the function exp(-iuit) . 
The electromagnetic strength E and B can be written as follows, 

E = E(r)exp(-iut) , B = B(r)exp(-iu;t) , (37) 

substituting the above in Eq.(36), 

{-(w/t;o) 2 + Sa|}E(r) = , {-(co/v ) 2 + Ed]} B(r) = . (38) 

From the Maxwell equations, <9 B + V* x E/V = and d E/Vb + V* x B = 0, we find that the 
amplitude of electromagnetic wave will be increased steadily or decreased continuously, although 
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this variation may be quite slow. So the field strength is the hyperbolic cosine cos(ia), and can be 
replaced by exp(-ila). 

And then, the field strength E(r) and B(r) are 

E(r) = E o exp(-ila) , B(r) = B o exp(-ila) , (39) 

where Eo and Bo both are the constant vectors in the octonion space. 
Further, substituting Eq.(39) in Eq.(38), we have the result, 

- (lu/v ) 2 + XK] = . (40) 

In the electromagnetic field, we will find that the electromagnetic waves are the transverse waves 
in a vacuum. From the field equations 

V • E = , V ■ B = , (41) 

we obtain 

K • E = , K • E' = , K • B = , K • B' = , (42) 

where E = Eo o I , B = Bq o I . The wave vector K = Y>(IjKj) . 

The above states that the electromagnetic waves, Eo and Bo , belong to the transverse waves. 
So do the two new wave components, E and B . Moreover, the amplitudes of wave components, 
E' and B , are equal to that of Eo and Bo respectively. 

In the electromagnetic field, there exist some relationships among the wave components, Eo and 
Bo , with two new kinds of wave components, E and B' . The field equations, 

VxE = k v dB/dt , V-E = , (43) 

can be rewritten as, 

V o E = k v dB/dt . (44) 

where the coefficient k v = Vo/vo, with k v ps 1 in general. 
Expanding the exp(-ila), the above implies that 

[K o E + cok v B ]cosh(a) + [K o E - ujk v B' ]sinh(a) = , (45) 

considering Eq.(42) and the a may be any value, and then 

K x E + k v uB = , K x E - k v uB' = . (46) 

In the same way, from the field equations 

voVqV x B = dE/dt , V • B = , (47) 

we have 

K x B + w(E /V r o)/« = , K x B - u(E' /V )/v = . (48) 

In the above, the K and Eo will yield a new wave component B' , while the K and Bo produce 
a new wave component E . Associating with the K x Bo is only the E , and with the K x Eo is 
the B' . This implies that it is impossible to determine simultaneously the Eo and Bo via the K. 
However, we can measure the Eo and B synchronously, or the E and Bo at the same time. There 
are many electric substances being able to slow down the speed of electromagnetic waves, and the 
field strength will possess the wave features, including the reflection and refraction etc. 



8 



3.4. Reflection and refraction 

In the electromagnetic field, by means of above features and the boundary condition, we can achieve 
the law of reflection, law of refraction, Fresnel formula, and the total internal reflection, etc. 

The electromagnetic wave transmits from the first medium to the second one, while there is one 
interface between these two media. In this interface, the t is the vector unit along to a tangent 
direction. The tangent components E< of E, and B[ of B' will be continuous respectively, for there 
is not electric current in the interface. There are the incident wave B" and reflective wave B^ in 
the first medium, and the refractive wave Bg in the second medium, with H' = B'/^ • Therefore 
the boundary condition in the interface are 

E a t + E\ = Et , [Bf + B*)/^ = B' t c /» 2 , (49) 

where E;f = tE( , Bj = tB' t , with f = a,b,c . The / = a stands for the incident wave, meanwhile 
f = b and / = c for the reflective and refractive waves respectively. [i\ and E\ are the permeability 
and permittivity for the first medium respectively. While [i 2 and e 2 for the second medium. 

In the Cartesian coordinate system oxyz, setting the plane z = on the interface, and the 
incident wave on the plane y = , the z-axis points to the second medium. In the incident plane, 
there is y = , and then K y = . In the interface, there exists z = 0, therefore 

E% t exp{ilxK^)exp{iuj a t) + E\ t exp\i\{xK b x + yK b )]exp(iu b t) 
E c ot exp\il{xK c x + yK c y )}exp{iLo c t) , (50) 
B' a t exp(axK^)exp(iuj a t)/ 'm + B$ t exp\iI(xK b + yK b y )]exp{iuj b t) / 'm 
= B' c t exp[iI(xK c x + yK c y )]exp{iu c t)/^ 2 , (51) 

where the wave vector K-^ = IiKl + I 2 Ky + I3.fr/ , the wave number K$ = cu? /Vq . 
Equating exponential parts for each coordinate and time, 

xK a x = xK b x + yK b y = xK% + yK c y , co a t = u b t = u c t , (52) 

and then 

K% = K b = K% , 0J a = J = u c , K b = K c y . (53) 

The above means that the angular frequencies of reflective wave and of refractive wave both are 
equal to that of incident wave. The refracted ray and reflected ray both are in the plane y = , 
that is these two rays and incident ray stand in same one plane. In the result, we will neglect the 
subscript (t) in the following context. 

According to the above results and definition of refractive index, we get 

K x = K a sin9 a , K b = K b sin9 b , K c x = K c sin9 c , (54) 

therefore 

6 a = 6 b , sin9 a /sin9 c = n 21 , (55) 

where 6 a , 9 b , and 9 C are the angle of incidence, angle of reflection, and angle of refraction respec- 
tively. n 2 \ is the refractivity of the second medium relative to the first one. 

The above results are so-called law of reflection and law of refraction for the light. 

3.5. Fresnel formula 

In the electromagnetic field, Fresnel formula can be derived from the relationships among the 
amplitudes and angles of incident wave, reflective wave, and refractive wave. 
The boundary condition of two media interface can be reduced to 

ES + E b = E c , (Bq + B' b )/^ = B' c /v 2 . (56) 

For the high-frequency vibration of electromagnetic waves, the fi of each media is almost equal 
to jjlq . And then the refractive index is 



n 2 \ = sin9 a /sin9 c = \fe 2 j~E\ 



(57) 
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When the E a is perpendicular to the incident plane, Eq.(46) says that 

B'o = ^e 1 ii Eq , Bq = ^/eiHoEq , B'q = ^e 2 ^Eq , (58) 
combining this with Eq.(56), gives 

E% + E% = E%, - y/e-[E%cos6 a + ^E b cos9 b = -^E c Q cos6 c , (59) 

and then 

E%/E$ = -sin(9 a - 9 c )/sin(9 a + 9 C ) , E%/E% = 2cos9 a sin9 c / sin{9 a + 6 C ) . (60) 

In the case of the E a is parallel to the incident plane, we have 

E%cos9 a - E b cos9 b = E c cos9 c , y/eiE$ + y/e~{E b = y/e5E% , (61) 

substituting Eq.(57) in the above, therefore 

E%/E% = tg(9 a - 9 c )/tg{9 a + 9 C ) , E%/E% = 2cos9 a sin9 c / sin(9 a + 9 c )cos(9 a + 9 C ) . (62) 

In the electromagnetic field, from Eqs.(46), (49), and (57), we obtain Fresnel formula, Eqs.(60) 
and (62), related to the E . On the other hand, we can find similar inferences about the E' , and 
some other results including Brewster's angle and the half-wave loss etc. 

3.6. Total internal reflection 

When the light passes from a denser medium to the interface of the less dense medium, e± > £2 , 
n2i < 1 , and 9 C > 9 a . Increasing the incident angle 9 a to the critical angle, the 9 C will achieve 
7r/2 . Meanwhile the refracted ray sweeps along the interface, and n 2 \ = sin9 a = \Jeils\ . 
Enlarging the incident angle continuously, we will find that sin9 a > n 2 \ , and then 

sin9 c = sin9 a /n 21 > 1 , cos9 c = i^(sin9 c ) 2 - 1 , (63) 

that means the refractive angle 9 C becomes an imaginary number. At the moment, the x-component 
and z-component of refracted wave vector K c are respectively, 

K% = K c sin9 c > K c , K c z = ±iy/ ' {K c x ) 2 - (K c ) 2 . (64) 
Considering K c £ = , and H(Kjrj) = xK® + zK c z , the electromagnetic strength will be 

E c = Eg o exp[-zI^(K^) 2 - {K c ) 2 ] o exp{-ilxK*)exp{iujt) , (65) 
B c = B c o exp[-zIy/(K£) 2 - (K c ) 2 } o exp(-iIxK^)exp{iujt) . (66) 

It is easy to find that the exp[—zl^/ (K£) 2 — (K c ) 2 ] is an octonion quantity. The above implies 
that the electromagnetic transmitting along the z-axis will be attenuated quickly. Eventually the 
incident wave is totally reflected inside. And that the amplitude and energy flux density of reflective 
wave are equal to that of incident wave respectively. 



Table 6: Some wave equations in the electromagnetic and gravitational fields. 

definition meaning condition wave equation separability 

X field quantity 

A = o X field potential A = 0* o (0 o X) = no 

B = o A field strength B = 0* o (0 o A) = no 

fiS = -{B/vq + 0)* oB field source (jB + B* o M/v Q = o (0* o B) = yes 

P = fiS/fMg linear momentum 

L = (K + k rx l€) o P angular momentum 

W = wo(B/wo + 0)oL energy and torque W-BoL = 0*o(0oL)=0 no 

N = w (B/u + 0)*oW power and force N - B* o W = ° (0* ° W) = yes 

C = v (B/v Q + 0)* o N physical quantity C - B* o N = o (0* o N) = yes 
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3.7. Other wave equations 

In the octonion space, the field potential A is defined from the field quantity X, which can be 
considered as the integral of field potential. And we can define the angular momentum L, energy- 
torque quantity W, power-force quantity N, and physical quantity C etc, from the radius vector M, 
field quantity X, and linear momentum P. 

Similarly there are some other kinds of wave equations in the electromagnetic and gravitational 
fields, including that for the field potential and angular momentum etc. 

4. TORQUE WAVE AND FORCE WAVE IN GRAVITATIONAL FIELD 

To incorporate various kinds of energies within a single definition, the angular momentum and 
energy will both be extended to apply within gravitational fields. The angular momentum density 
= ^0 + is defined from the linear momentum density ¥ g = fi g Eg/n g in Eq.(l) and the 

radius vector R 9 . And it can be rewritten as follows, 

L 9 = (Rg + k rx X g ) o Fg , (67) 

where Z = (r + k rx x )po + (r + k rx x) ■ p , 1 = (r + k rx x )p + p (r + k rx x) + (r + k rx x) x p, and 
1 = T,(ljij). The quaternion quantity X ff = E(ccjij) is similar to Hertz vector in the electrodynamics 
theory, with k rx X g <C M ff and k rx = 1 . The derivation of X g will yield the gravitational potential. 
The do = vq(8oXo + V • x) and a = t>o(<9o x + Vxq + V x x) are the scalar potential and vectorial 
potential of gravitational field respectively, x = T,(xjij). 

We choose the following definition of energy to include various energies in the gravitational field. 
In quaternion space, the quaternion energy-torque density W g = wo + H(wjij) is defined from the 
angular momentum density h g , 

W g = v (M g /v + 0)ohg , (68) 

where u>o = fo^Mo + ^oV • 1 + (g/^o + b) ■ 1, and the —wo/2 is the energy density, including the 
density of kinetic energy, potential energy, field energy, and work etc. While the w = Y*(wjij), with 
w = vodol + voVlo + t>oV x 1 + lo(g/vo + b) + (g/^o + b) x 1. And the — w/2 is the torque density, 
including that caused by gravity and some other force terms. 

The quaternion power- force density N 9 = no + T,(rijij) is defined from the W g , 

Ng = v (B 9 /v + Oy oWg , (69) 

where no = vodowo + t>oV* ■ w + (g/^o + b)* • w, and /q = —no/(2vo) is the power density. While 
n = E(rijij) = voV*wo + -uo^ow + ^oV* x w + (g/t>o + b)* x w + wo(g/v + b)*. The force density 
f in the gravitational field can be defined from n , that is, f = — n/(2uo) • The force density f 
includes the density of gravity, inertial force, and energy gradient etc. 

To deduce the wave equation for force, the physical quantity C 9 = Co + T,(cjij) can similarly be 
defined as follows, 

Cg=V (Mg/v + 0)* °Ng , (70) 

where cq = v^don^ + vqV* ■ n + (g/t>o + b)* • n, while c = S(cjij) = vqV*uq + vq8qii + vqV* x n + 
(s/ v o + b)* x n + no(g/vQ + b)*. The c includes the derivative of acceleration (or jerk). 

4.1. Torque wave equation 

Choosing the gauge equation, wo = 0, from Eq.(68) to simplify following calculation. The torque 
density w can be separated into three components, w = w e /wo + w;, + w s , 

w e /v = v d l + vqVIo , w b = v V x 1 , (71) 
w 8 = i (g/wo + b) + (g/«o + b)xl. (72) 

4.1.1. Case A 

Similar to Eqs.(2) and (3), in case of wq = and w s = 0, the definition equation of quaternion 
power-torque density, Eq.(69), can be decomposed into the Maxwell-like equations, 



V • (w 6 /«o) = , d (w b /v ) + V* x (w e /v )/v = , 
V* • (w e /vo)/vo = uq/vq , d (-w e /v ) /v + V* x (-Wb/vo) = n/v 



(73) 
(74) 
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In the above, the w s = is equivalent to g/vo + b = from Eq.(72) when L^O. If there exists 
further N 9 = , the above or Eq.(69) can be reduced to, 



0* O Wg = 

Applying the operator to the above, we get, 



. 



(d 2 + d\ + dl + d! 



. 



(75) 



(76) 



The above is Laplace equation as well as the wave equation for the torque in the quaternion 
space. Proceeding with the operator do and V, we can obtain the wave equation about the torque 
components from Eqs.(73) and (74) directly, 



{dl + d\ + dl + dl)w e = o , {dl + dl + d\ + <9f)w& = o . 



(77) 



The above means that two components of torque, w e and w&, are both possessed of the wave 
features in the gravitational field. And these torque component is similar to the periodic torque in 
the vibrational mechanics partially, and may cause the vortex street in the fluid mechanics. 



Table 7: Maxwell-like equations for field potential and for angular momentum. 



field potential 



angular momentum 



field strength 
gauge equation 
field equation 



g/v = d a + Va 
b = V x a 
&o = 



-fjSg = (Vg/VO + 0)* 0Mg 

V = 

Maxwell-like equations V • b = 

d b + V* x g/v Q = 
V* • g = -ml Eg 
dpg/^o + V* x b = -/j, g s 



w e /vo = v d l + v Vl 
w;, = v V x 1 
w = 
w s = 

-LlNg = {Mg/V + 0)* O Wg 

/i = -l/v 

V • (w 6 /«o) = 

<9 (w b /v ) + V* x (w e /vo)/v a = 
V* • {w e /v )/vo = rio/vo 
d {w e /v )/vo + V* x (w 6 /t> ) = n/u 



4.1.2. Case B 

In case of wq = and w s / 0, Eq.(69) can be decomposed as follows, 

n /v = V* • w + (g/uo + b)* • w/v , (78) 
ii/vq = d w + V* x w + {g/v + b)* x w/v . (79) 

In general the above can be separated into the equations more complex than Eqs.(73) and (74), 
and may possess more complicated wave features than which in Eqs.(76) and (77). 
Specially there may exist following conditions, 

V* • (w s + w 6 ) + {g/v + b)* • w/v = , (80) 
3 (w s + w 6 ) + V* x (w e /v + w 8 ) + {g/v + b)* x w/v = , (81) 

and considering Eq.(73), the above conditions can be reduced to, 

V* • w s + {g/v + b)* • w/v = , (82) 
a w s + V* x w s + {g/v + b)* x w/v = . (83) 

Under the above conditions when w s / 0, Eqs.(78) and (79) can be reduced to the Maxwell-like 
equations, which is same as Eqs.(73) and (74). 

If there exists additionally N 9 = , Eqs.(78) and (79) can further be simplified to, 

0* o W 9 = , (84) 

and then we have Eqs.(76) and (77), and some similar inferences about the torque components. 

In Eq. (69) , when there exists the condition relation N g — M* o W s = , we should gain some 
more complex deductions than the above directly. 
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4-1.3. Transverse wave 

For the torque wave with angular frequency u, the torque wave should be a harmonic function 
cosut , and then should be chosen as the function exp(-iujt) . The torque wave w e and w b can be 
written as follows, 

w e = w e (r)exp(— iuit) , w& = w b (r)exp(— loot) , (85) 
combining this with Eq.(76), gives 

{-( w /t, ) 2 + Sa|}w e (r) = , {-(w/t;o) 2 + Sa|}w5(r) = . (86) 

From the analysis of Eqs.(73) and (74), we find that the torque wave should be one hyperbolic 
cosine cos{if3) , and then should be chosen as the function exp(—i\(3) . The wavedike vector 
k w = Y>(ijk w j) , while the j3 = —Y,(rjk w j) , with k w j being the coefficient. 

The torque wave w e (r) and wj(r) are 

w e (r) = w e0 o exp(-ii/3) , w b (r) = w 50 o exp(-ii/3) , (87) 

where w e o and w^o both are constant vectors in the quaternion space. 
And then we have the result from the above and Eq.(85) , 

- (oj/vo) 2 + Xk 2 wj = . (88) 

In the gravitational field, we will find that torque waves are the transverse waves in a vacuum. 
Similar to Eqs.(16)~(23), there are following results from Eqs.(73) and (74), 

k w • w e0 = , k w ■ Wg = , k w ■ w b0 = , k w ■ w' 60 = , (89) 
k w x w' e0 + ww fe0 = , k w x w e0 - ujW bQ = , (90) 
k w x w^ + uj(w e0 /v )/v = , k w x w b0 - uj(w' e0 /v )/v = , (91) 

where w e0 = w e0 o i , w^ = w b0 o i . 

The above means that there are some relationships between W 9 with W g o i . The k w and w e o 
will yield a new component w' M , while the k w and W50 produce the wj, . All of these torque 
components belong to the transverse waves. Moreover, the amplitudes of torque wave components, 
Wg and w^ , are equal to that of w e o and w^o respectively. 

4.2. Force wave equation 

In a similar way to the case of torque wave, we can study some wave features about the force in the 
gravitational field. In the quaternion space, we choose the mass continuity equation as the gauge 
equation from Eq.(69) to simplify the following calculation, that is uq = 0. 

The force density n can be separated into three components, n = n e /vo + + n s , 

n e /v = v d w + v Vw , n 6 = v V x w , (92) 
n s = w {g/v + b) + (g/v + b) x w . (93) 

4.2.1. Case A 

Similar to Eqs.(73) and (74), in case of no = and n s = 0, the definition equation of quaternion 
physical quantity, Eq.(70), can be decomposed into the Maxwelldike equations, 

V • (n b /v ) = , d {n b /v ) + V* x {n e /v )/v = , (94) 
V* • (n e /v )/v = c /v , d (n e /v )/v + V* x (n b /v ) = c/v . (95) 

In the above, the n s = is equivalent to g/vQ + b = from Eq.(93) when W 9 / 0. If there 
exists further B g = and C g = , the above or Eq.(70) can be reduced to, 

0* o n g = . (96) 

Applying the operator to the above, we get the wave equation 

(dl + d\ + d 2 2 + dl)n g = . (97) 



13 



The above is Laplace equation for the force in the quaternion space too. Proceeding with the 
operator do and V, we can obtain the wave equation about the force components, n e and n b , from 
Eqs.(94) and (95), 

{dl + df + d\ + df K = o , {dl + d\ + d\ + df H = o . (98) 

The above means that two components of force, n e and n&, are both possessed of the wave 
features in the gravitational field. These force component is similar to the periodic force in the 
vibrational mechanics partially, and may activize the vortex street in the fluid mechanics. 



Table 8: Maxwcll-likc equations for field potential and for torque. 



field potential 



torque 



field strength 
gauge equation 
field equation 



g/w = <9 a + Va 
b = V x a 
b = 



-/jS g = (Vd + 0)*oB 9 
V = Mg 

Maxwell-like equations V • b = 

d h + V* x g/v = 
V* • g = - fn/sg 

dog/v + V* x b = -ii g s 



n e /v = v dow + voVw 
rib = «oV x w 
n = 
n s = 

-/iC 9 = (B 9 /« + 0)*oN 9 

A* = -1M) 

V • (n b /v ) = 

d (n b /v ) + V* x (n e /v )/v = 
V* • (n e /v )/v = c /v 
do(n e /vo)/v + V* x (n b /v ) = c/v 



4.2.2. Case B 

In case of no = and n s / 0, Eq.(70) can be decomposed as follows, 

c /v = V* • n + (g/uo + b)* • n/v , (99) 
c/v = d n + V* x n + (g/v + b)* x n/v . (100) 

In general the above can be separated into the equations more complex than Eqs.(94) and (95), 
and may possess more complicated wave features than which in Eqs.(97) and (98). 
Specially there may exist following conditions, 

V* • (n s + n b ) + (g/«o + b)* • n/v = , (101) 
d (n s + n b ) + V* x {n e /v + n 8 ) + {g/v + b)* x n/v = , (102) 

and considering Eq.(94), the above conditions can be reduced to, 

V* • n s + (g/v + b)* • n/v = , (103) 
d n s + V* x n s + (g/v + b)* x n/v = . (104) 

Under the above conditions when n s / 0, Eqs.(99) and (100) can be simplified to the Maxwell- 
like equations, which is same as Eqs.(94) and (95). 

If there exists additionally C 9 = , Eqs.(99) and (100) can further be reduced to, 

0* o N 9 = , (105) 

and then we have Eqs.(97) and (98), and some similar inferences about the force components. 

In Eq.(70), when there exists the condition relation C g — B* o N 9 = , we should gain some 
deductions more complex than the above directly. And the alternating field strength, torque waves, 
and force waves etc may cause jointly the density waves in the galactic dynamics of astronomy. 

4-2.3. Transverse wave 

For the force wave with angular frequency cu, the force wave should be a harmonic function cosivt 
as well, and then should be chosen as the function exp{—iuit) . The force wave n e and can be 
written as follows, 



n e = n e (r)exp(—iujt) , n^, = nb(r)exp{—iojt) , 



(106) 
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combining this with Eq.(97), gives 

{-(cj/v ) 2 + i:d]}n e (r) = , {-(uj/v ) 2 + Zd]}n b (r) = . (107) 

From the analysis of Eqs.(94) and (95), we find that the force wave should be one hyperbolic 
cosine cos(ij) as well, and then should be chosen as the function exp(-ii^f) . The wave-like vector 
k n = Y>(ijk n j) , while the 7 = —Y>(rjk n j) , with k n j being the coefficient. 

The force wave n e (r) and n&(r) are 

n e (r) = n e0 o exp(-ii-y) , n b (r) = n b0 o exp(-ii-y) , (108) 

where n e o and n^o both are constant vectors in the quaternion space. 
And then we have the result from the above and Eq.(106) , 

- (oj/v ) 2 + £A£ . = o _ (109) 

In the gravitational field, we will find that force waves are the transverse waves in a vacuum. 
Similar to Eqs.(89)~(91), there are following results from Eqs.(94) and (95), 

k n • n e0 = , k n • n^ = , k n • n 60 = , k n • n' M = , (110) 
k n x n' e0 + ujn b0 = , k n x n e0 - un' b0 = , (111) 
k„ x n^ + u(n eQ /vo)/v = , k n x n 60 - Lu(n' e0 /vo)/v = , (112) 

where n' e0 = n e0 o i , n' 60 = n 60 o i . 

The above means that there are some relationships between N 9 with N 9 o i . The k ra and n e o will 
yield a new component n^ , while the k n and n^o produce the n' e0 . All of these force components 
belong to the transverse waves. Moreover, the amplitudes of force wave components, n' e0 and n^ , 
are the same as that of n e o and n^o respectively. 



Table 9: The physical quantity and related waves in the gravitational field. 





field strength 


torque 


force 


static quantity 


yes 


yes 


yes 


alternating quantity 


yes 


yes 


yes 


independent wave 


yes 


yes 


yes 



5. TORQUE WAVE AND FORCE WAVE IN ELECTROMAGNETIC FIELD 

The torque wave feature as well as the force wave feature can be described simultaneously by the 
algebra of octonions in the gravitational field and electromagnetic field. 

In the case for coexistence of electromagnetic field and gravitational field, the octonion angular 
momentum density L = S(Zjtj + Lj/j) is defined from the octonion radius vector R and octonion 
field source S in Eq.(24). 

L = (R + k rx X) o P , (113) 

with 

k = (ro + k rx x )po + (Ro + k r xX ) o P + (r + k rx x) • p + (R + k rx X) ■ P , 
1 = (r + k rx x )p + (R + k rx X) o P + (r + /c ra x) x p + (R + k rx X ) o P 
+p (r + k rx x) + (R + k rx X) x P , 
Lo = (?"o + k rx x )~Po + (r + k rx x) ■ P + p (Ro + k rx X ) + (R + k rx X) ■ p , 
L = (r + k rx x )~P + (r + k rx x) o P + (r + k rx x) x P + (R + k rx X ) o p 
+Po(R + k rx X) + (R + k rx K) x p , 

where P = fjS/n g . L = L I , L = Z(Ljlj). P = P I , P = Z(Pjlj). R = R I , R = ^{Rjlj). 
X = Yi(xiii + Xili). Xo = ^o-^O) ^ = ^(Xjlj). The derivation of octonion physical quantity X 
yields the gravitational potential as well as the electromagnetic potential. 
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The octonion energy-torque density W = + WjJj) can be denned from the octonion field 

strength B and octonion angular momentum density L in Eq.(113). 

W = ti (BM + 0)oL, (114) 

with 

w = v d l + v V ■ 1 + (g/v + b) • 1 + k eg (E/V + B) ■ L , 
w = v d l + v X7l + v V x 1 + l Q (g/v Q + b) + (g/v + b) x 1 
+k eg (E/V + B) x L + k eg (E/V + B) o L , 
W = v d U + k eg (E/V + B) • 1 + v V ■ L + (g/v + b) • L , 
W = v d L + l k eg {E/V + B) + v V x L + v V o L + k eg (E/V + B) x 1 
+(g/v + b) x L + (g/«o + b) o L , 

where —wo/2 is the energy density, includes the density of kinetic energy, potential energy, field 
energy, work, electric potential energy, magnetic potential energy, and interacting energy between 
dipole moment with electromagnetic strength etc. w/2 = T,(wjij)/2 is the torque density, includes 
that caused by gravitational force, electromagnetic force, and other force terms etc. Wo = Woibj 
W = T,{WjIj). A + A = k eg A e . A = WoXo + V • X) and A = V (d X + V o X + V x X) are 
related to the scalar potential and vectorial potential of electromagnetic field respectively. 

The octonion power-force density N = £(njij + iVjJj) is defined from the octonion energy-torque 
density W in Eq.(114) and field strength B , 

N = vo(B/wo + 0)*oW , (115) 

with 

n = v d w + v V* • w + (g/uo + b)* • w + k eg (E/V + B)* • W , 
n = v d w + v V*w + v V* x w + ifo(g/^o + b)* + (g/v + b)* x w 
+k eg (E/V + B)* x W + k eg (E/V + B)* o W , 
N = v d W + k eg (E/V + B)* • w + t; V* • W + (g/v + b)* • W , 
N = v d W + w k eg (E/V + B)* + v V* x W + v V* o W 

+k eg (E/V + B)* x w + (g/« + b)* x W + (g/uo + b)* o W , 

where No = NqIq, N = Y>(NjIj). The /o = — no/(2vo) is the power density, while the vectorial part 
n = T,(rijij) is the function of forces. The f = — n/(2vo) is the force density, which includes the 
density of inertial force, gravity, Lorentz force, gradient of energy, and interacting force between 
dipole moment with magnetic strength etc. 

The octonion physical quantity C = S(cjtj + CjJj) is defined from the octonion power- force 
density N in Eq.(115) and field strength B , 

C = v (M/v + 0)* oN , (116) 

with 

c = v d n + v V* ■ n + (g/v + b)* • n + k eg (E/V + B)* • N , 
c = t>o<9on + v V*n + v V* x n + n (g/v + b)* + (g/v + b)* x n 
+k eg (E/V + B)* x N + k eg (E/V + B)* o N , 
C = v d N + k eg (E/V + B)* • n + t; V* • N + {g/v + b)* • N , 
C = v d U + n k eg (E/V + B)* + v V* x N + t; V* o N 

+k eg (E/V + B)* x n + (g/v IbfxNl (g/v + b)* o N , 

where Co = Colo, C = T,(CjIj). 
5.1. Torque wave equation 

Comparing with the classical electromagnetic theory, there are the W ff and N 9 in the E 9 space, 
while the W e and N e in the E e space. And the physical quantity W as well as N are devised to 
describe consistently in the electromagnetic and gravitational fields. 
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According to the basis vectors, Eq.(114) can be decomposed further as follows, 

(B; o W g + k eg W e o W e )/v + 0* o Wg = N g /vo , (117) 

(B* o W e + k eg M* e o W g )/V + 0* O W e = Ne/^0 , (118) 

where W e = ^{WJi) , N e = S(JVjJj) . In the above, Eq.(117) is similar to Eq.(69) and is suitable 
for the torque wave equation, while Eq.(118) is for the torque-like wave equation. 

Choosing the gauge equation, wq = 0, from Eq.(114) to simplify following calculation. The 
torque density w can be separated into three components, w = w e /vo + Wf, + w s , 

w e /v = v d l + v Vl , w b = v V x 1 , (119) 
w s = Z (g/u + b) + {g/v + b) x 1 + k eg {B/V + B) x L + k eg {B/V + B) o L . (120) 

Substituting Eqs.(117), (119), and (120) for Eqs.(69), (71), and (72) respectively, there are some 
deductions about the torque wave features similar to Eqs.(73)~(91) in the case for coexistence of 
gravitational and electromagnetic fields. 

5.2. Torque-like wave equation 

From Eq.(114) in the octonion space, we choose the gauge equation, Wo = 0. And the torque-like 
quantity W can be separated into three components, W = W e /Vo + Wj + W s , 

We/Vb = v d L + ^ V o L , W 6 = v V x L , (121) 
W s = l k eg {B/V + B) + (g/v + b) o L + k eg (E/V + B) x 1 + {g/v + b) x L . (122) 

5.2.1. Case A 

Similar to Eqs.(73) and (74), in case of Wo = and W s = 0, Eq.(118) can be decomposed into 
the Maxwell-like equations, 

V • (W b /v ) = , d (W b /v ) + V* x (W e /V )/v = , (123) 
V* • (W e /V )/v = No/vo , d (W e /V )/v + V* x (W h /v ) = N/i> . (124) 

In the above, the W s = means that g/vQ + h = and E/Vo + B = from Eq.(122) when 
L / 0. If there exists N e = , the above or Eq.(118) can further be reduced to, 

0* o We = . (125) 

Applying the operator to the above, we have the wave equation or Laplace equation 

o (0* o W e ) = , or (dl + d\ + d\ + <9f )W e = . (126) 

Proceeding with the operator do and V, we can obtain the wave equation about the components 
of torque-like quantity from Eqs.(123) and (124), 

{dl + dl + d 2 2 + dl)w e = o , (d 2 + d\ + dl + d 2 3 )w b = o . (127) 

The above means that two components of torque- like quantity, W e and W5, are both possessed 
of wave features in the electromagnetic field. 

5.2.2. Case B 

In case of Wo = and W s 7^ 0, Eq.(118) can be decomposed as follows, 

No/t>o = V* • W + (g/«o + b)* • W/^o + k eg (E/V + B)* • w/v , (128) 
N/t>o = d W + V* x W + (g/uo + b)* x W/t; 

+k eg (E/V + B)* x w/v + w k eg (E/V + B)*/v . (129) 

In general the above can be separated into the equations more complex than Eqs.(123) and 
(124), and may possess more complicated wave features than which in Eqs.(126) and (127). 
Specially there may exist following conditions, 

0= V*-(W 8 + W 6 ) + (g/«o + b)*-W/«o + fce fl (E/yo + B)*-w/«o , (130) 
= d (W s + W 6 ) + V* x (We/Vb + W s ) + (g/v + b)* x W/t; 

+k eg {B/V + B)* x w/v + w k eg (E/V + B)7«o , (131) 



17 



Tabic 10: Maxwell-like equations for torque and torque-like in the electromagnetic field. 

torque torque-like 

field strength w e /w = i>o<9ol + v Vlo W e /Vb = v doL + w V o L 

w 6 = v \7 x 1 W fc = w V x L 

gauge equation w = W = 

w s = W s = 

field equation N 9 =i;oW 9 | /c eg B* o W e + w 0* o W 9 N e = 1* o W e + k eg M* e o W s + w 0* ° W e 

(J,= -1/V fJ, = -l/v 

Maxwell-like V • (w b /v ) = V • (W b /v ) = 
equations 5 (w 6 /w ) + V* x {w e /v )/v = <9o(W 6 /w ) + V* x {W e /V )/v = 
V* • {w e /v )/v = n /v V* • {W e /V )/v = N /v Q 
do(w e /vo)/v + V* x {wh/vp) = n/vp d {W e /V )/v + V* x {W b /v ) = N/v a 

and considering Eq.(123), the above conditions can be reduced to, 

0= V* • W s + (g/^o + b)* -W/vo + kegiE/Vo + B)* --w/vo , (132) 
= 8 W S + V* x W s + (g/uo + b)* x W/w + /e e9 (E/y + B)* x w/t;o 

+w A;e ff (E/^o + B)7«o , (133) 

Under the above conditions, from Eqs.(128) and (129) when W s / 0, we can obtain the Maxwell- 
like equations, which are same as Eqs.(123) and (124). 

If there exists additionally N e = , Eqs.(128) and (129) can further be reduced to, 

0* o W e = , (134) 

and then we have Eqs.(126) and (127), and some similar inferences about torque-like components. 

In Eq.(118), when there exists the condition relation N e - (B* o W e + k eg M* e o W 9 ) = , we 
should gain some more complex deductions than the above directly. 

5.2.3. Transverse wave 

For the torque wave with angular frequency u, the torque wave should be a harmonic function 
cosut , and then should be chosen as the function exp(—iojt) . The torque wave W e and W b can 
be written as follows, 

W e = W e {r)exp{-iojt) , W b = W b {r)exp{-iujt) , (135) 
combining this with Eq.(126), gives 

{-(^ ) 2 + Sa|}W e (r) = , {-(u/v ) 2 + Xd]}W b (r) = . (136) 

From the analysis of Eqs.(123) and (124), we find that the torque wave should be one hyperbolic 
cosine cos{ij3) , and be chosen as the function exp(— il/3) . The wave-like vector K. w = T,(IjK w j) , 
while the ft = —T,(rjK w j) , with K w j being the coefficient. 

The torque wave W e (r) and W&(r) are 

W e (r) = W e0 o exp(-ilp) , W 6 (r) = W 60 o exp(-ilp) , (137) 

where W e o and Wjj both are constant vectors in the octonion space. 
And then we have the result from the above and Eq.(135) , 

- (lj/vo) 2 + VKlj = . (138) 

In the electromagnetic field, we will find that torque waves are the transverse waves in a vacuum. 
Similar to Eqs.(110)~(112), there are following results from Eqs.(123) and (124), 

K,i, • W e0 = , K w ■ W' e0 = , K w ■ W b0 = , K w ■ W' 60 = , (139) 
K w x W' e0 + kyuWto = , K„, x W e0 - k v ujW' m = , (140) 
K w x W' b0 + co(W e0 /V )/v = , K w x W b0 - u(W' e0 /V )/v = , (141) 

where W e0 = W e0 o I , W' 60 = W 60 o I . 

The above means that there are some relationships between W e with W e o I . The JZ W and W e o 
will yield a new component W^, , while the K w and Wfto produce the W e0 . All of these torque 
components belong to the transverse waves. Moreover, the amplitudes of torque wave components, 
W' e0 and W£, , are the same as that of W e o and W^q respectively. 
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5.3. Force wave equation 

Similar to Eq.(115), the C is devised to describe the octonion physical quantity consistently in 
the electromagnetic field and gravitational field. According to the basis vectors, Eq.(116) can be 
decomposed further as follows, 

(i; o N 9 + k eg W e o N e )/v + 0* o N 9 = Cg/vo , (142) 
(i; o N e + k eg W e o n g )/v + 0* o N e = C e /v , (143) 

where C e = S(CjJj) . In the above, Eq.(142) is similar to Eq.(117) and is suitable for the force 
wave equation, while Eq.(143) is for the force-like wave equation. 

To simplify the following calculation, we choose the mass continuity equation from Eq.(115) as 
the gauge equation in the electromagnetic field, that is no = . The force density n can be divided 
into three components, n = n e /vo + + n s , 

n e /v = v d w + v X7*w , n b = v V* x w, (144) 
n s = wo(g/vo + b)* + (g/u + b)* x w + k eg {B/V + B)* x W + k eg {B/V + B)* o W .(145) 

Substituting Eqs.(142), (144), and (145) for Eqs.(70), (92), and (93) respectively, there are some 
deductions about the force wave features similar to Eqs.(94)~(112) in the case for coexistence of 
gravitational and electromagnetic fields. 

5.4. Force- like wave equation 

From Eq.(115) in the case for coexistence of electromagnetic and gravitational fields, we choose the 
charge continuity equation as the gauge equation, that is No = 0. And the force-like quantity N 
can be separated into three components, N = N e /Vb + N& + N s , 

N e /V = v d W + v V* o W , N b = v V* x W, (146) 
N s = w k eg (E/V + B)* + (g/«o + b)* o W + k eg (E/V + B)* x w + (g/v + b)* x W.(147) 



Tabic 11: Maxwell- like equations for force and force- like in the electromagnetic field. 

force force-like 

field strength n e /v = v do\v + v Vwo N e /Vb = v doW + w V o W 

lib = woV x w N(, = u V x W 

gauge equation no = No = 

n s = N s = 

field equation C 9 = B* o N 9 + k eg M* e o N e + v <>* o N 9 C e = B* o N e + k eg M* e o N s + w 0* o N e 

H = -l/v n = -l/v 

Maxwell-like V • (n b /v ) = V • (N b /v Q ) = 

equations d (n b /v ) + V* x (n e /v a )/v = d (N h /v ) + V* x (N e /V )/v a = 

V* • (n e /v )/v Q = c a /v a V* • (N e /V )/v = C /v Q 
d (n e /v )/v Q + V* x (n b /v Q ) = c/v d (N e /V )/v Q + V* x (N b /v ) = C/vp 



5.4.I. Case A 

Similar to Eqs.(94) and (95), in case of No = and N s = 0, Eq.(143) can be decomposed into the 
Maxwell- like equations, 

V • (N 6 /«d) = , d (N b /v ) + V* x {N e /V )/v = , (148) 
V* • {N e /V )/v = Co/vo , d (N e /V )/v + V* x (N b /v ) = C/v . (149) 

In the above, the N s = means that g/v Q + b = and E/Vb + B = from Eq.(147) when 
W^O. If there exists C e = , the above or Eq.(143) can further be reduced to, 

0* o N e = . (150) 

Applying the operator to the above, we have the wave equation or Laplace equation 

o (0* o N e ) = , or {dl + d\ + d\ + a|)N e = . (151) 
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Proceeding with the operator do and V, we can obtain the wave equation about the components 
of force-like quantity from Eqs.(148) and (149), 

{8l + d\ + d\ + df )N e = , (dl + d\ + 8l + <9f )N 6 = . (152) 

The above means that two components of force-like quantity, N e and N;,, are both possessed of 
wave features in the electromagnetic field. 

54.2. Case B 

In case of No = and N s / 0, Eq.(143) can be decomposed as follows, 

C /v = V*-N+(g/t;o + b)*-N/t;o + A :e9 (E/yo + B)*-n/«o , (153) 
C/v = <9 N + V* x N + (g/«o + b)* x N/v 

+k eg (E/V + B)* x n/vo + n k eg (E/V + B)*/v . (154) 

In general the above can be separated into more complex equations than Eqs.(148) and (149), 
and may possess more complicated wave features than that in Eqs.(151) and (152). 
Specially there may exist following conditions, 

0= V* ■(N s + -N b ) + (g/v + h)* -N/vo + hgiE/Vo + B)* -n/v , (155) 
= 8 (N S + N b ) + V* x (NJV + N s ) + (g/v + b)* x N/v 

+k eg {E/V + B)* x n/v + n k eg {E/V + B)7«o , (156) 

and considering Eq.(148), the above conditions can be reduced to, 

= V* • N s + (g/v + b)* • N/v + k eg (E/V + B)* • n/v , (157) 
= 8 N S + V* x N s + (g/«o + b)* x N/t; + k eg {E/V + B)* x n/v 

+n k eg (E/V + B)*/v , (158) 

Under the above conditions, from Eqs.(153) and (154) when N s 7^ 0, we can obtain the Maxwell- 
like equations, which are same as Eqs.(148) and (149). 

If there exists additionally C e = , Eqs.(153) and (154) can further be reduced to, 

0* o N e = , (159) 

and then we have Eqs.(151) and (152), and some similar inferences about force-like components. 

In Eq.(143), when there exists the condition relation C e — (B* o N e + k eg W e o N 9 ) = , we should 
gain some more complex deductions than the above directly. 

In the case for coexistence of gravitational and electromagnetic fields, the alternating field 
strength, torque waves, torque-like waves, force waves, and force-like waves etc may agitate jointly 
the density waves in the galactic dynamics of astronomy. 



Tabic 12: The physical quantity and their waves in the electromagnetic field. 





field strength 


torque 


torque- like 


force 


force- like 


static quantity 


yes 


yes 


yes 


yes 


yes 


alternating quantity 


yes 


yes 


yes 


yes 


yes 


independent wave 


yes 


yes 


yes 


yes 


yes 



5.4-3. Transverse wave 

For the force wave with angular frequency u, the force wave should be a harmonic function cosujt , 
and then should be chosen as the function exp(-iujt) . The force wave N e and N5 can be written 
as follows, 

N e = N e (r)exp(-iujt) , N 6 = N b (r)exp(-iujt) , (160) 
combining this with Eq.(151), gives 

{-(w/t;o) 2 + Sa|}N e (r) = , {-(uj/v ) 2 + ^d]}N b (r) = . (161) 
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From the analysis of Eqs.(148) and (149), we find that the force wave should be one hyperbolic 
cosine cos{i~f) , and be chosen as the function exp(-il^) . The wave-like vector K n = Y>(IjK n j) , 
while the 7 = —T,(rjK n j) , with K n j being the coefficient. 

The force wave N e (r) and N&(r) are 

N e (r) = N e0 o exp(-il'y) , N 6 (r) = N 60 o exp(-ilj) , 

where N e o and N^q both are constant vectors in the octonion space. 
And then we have the result from the above and Eq.(160) , 

- (u;/v ) 2 + Y,Kl 3 = . (162) 

In the electromagnetic field, we will find that force waves are the transverse waves in a vacuum. 
Similar to Eqs.(110)~(112), there are following results from Eqs.(148) and (149), 

K n • N e0 = , K n ■ N' e0 = , K n ■ N b0 = , K n ■ N' M = , (163) 
K n x N' e0 + k v ujN b0 = , K n x N e0 - k v ooN' bo = , (164) 
K n x N' 60 + co(N e0 /V )/v = , K n x N 60 - u(N' e0 /V )/v = , (165) 

where N' e0 = N e0 o I , N' 60 = N 60 o I . 

The above means that there are some relationships between N e with N e o I . The K n and N e o 
will yield a new component N^ , while the K n and N^q produce the N' e0 . All of these force 
components belong to the transverse waves. Moreover, the amplitudes of force wave components, 
Ng and N' w , are the same as that of N e o and N^o respectively. 

6. FLUX WAVES IN ELECTROMAGNETIC AND GRAVITATIONAL FIELDS 

From the property of octonion spaces, we can define the linear momentum flux L p , angular momen- 
tum flux L/, energy-torque flux L w , power- force flux L n , and physical quantity flux L c respectively 
from the velocity V, field potential A, and linear momentum P etc. 

Consequently there may exist some other kinds of wave equations in the electromagnetic and 
gravitational fields for those flux quantities. 



Tabic 13: Energy flux an others in the electromagnetic and gravitational fields. 
linear momentum flux angular momentum flux energy — torque flux 
L p = (V + k rx A) o P L' = (V + k rx A) o L h w = (V + k rx A) o W 

W p = v (M/v + ())ohP W' =v {M/v + 0)o^ 1 W w =v (R/v + <>)oL w 
W = v {M/v + <»* oW p N' = v (M/v + 0)* o W ; r=«o(E/wo + 0)*oW'' 
CP = v (M/v a + 0)* oW C'=bo(B/bq + 0)'oW' C° = d (B/^ + 0)'oW" 



To incorporate various linear momentum flux within a single definition, the linear momentum 
flux L p is defined from the linear momentum density P = n g S/p, g in Eq.(24), the velocity V, and 
the field potential A . The related equations can be written as follows, 

LP = (V + k rx A) o P , W p = v (K g /v + 0) o Lf , (166) 
W = v {M g /v + 0)* o WP , CP = voIb 9 /v + 0)* o W . (167) 

From the above, we can describe the property of linear momentum flux in the electromagnetic 
field in the presence of gravitational field. 

We choose following definition of angular momentum flux to include various angular momentum 
flux in the electromagnetic field and gravitational field. In octonion space, the equations about the 
angular momentum flux are defined from the angular momentum density L , and the velocity V 
and field potential A, 

= (V + k rx A) oh, W l = v (M g /v + 0) o L' , (168) 
N z = v (M g /v + 0)* 0W 1 , C' = v (M g /v + 0)* o N z . (169) 

Similarly the octonion energy-torque flux equations of electromagnetic field and gravitational 
field are defined as follows, 

h w = (V + k rx A) o W , W w = v (M g /v + 0) o h w , (170) 
N w = v {Mg/v + 0)* o W w , C w = v {»g/v + 0)* o N w . (171) 
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The power-force flux and related equations in the octonion space can be defined from the power- 
force N , 

L n = (V + k rx A) o N , W n = v (M g /v + 0) o L™ , (172) 
N n = v (» g /v + 0)* o W n , C n = v (M g /v + 0)* o W l . (173) 

To deduce the wave equation for C flux, we can similarly define following equations, 

L c = (V + k rx A) o C , W c = v (K g /v + 0) o L c , (174) 
N c = v (M g /v + 0)* o W c , C c = v {Mg/v + 0)* o N c . (175) 

From the above we will deduce some wave equations and transverse waves about the flux waves 
similar to that of Eqs.(67)~(164), including the inferences about the energy flux. 

Table 14: Some flux quantities in the electromagnetic and gravitational fields. 
power — force flux C flux 

L" = (V + k rx A) o N L c = (V + k rx A) o C 

W" = v a (M/v a + 0) o L" W c = v a (M/v a + 0) o L c 
N™ = v (M/v + 0)* o W" N c = v (M/v + 0)* o W c 
C" = v (M/v Q + Q)* o N" C c = v Q (M/v + Q)* o N c 



Further, substituting (D + k rx M) for (V + k rx A) in Table 13 and Table 14, we will obtain some 
other sorts of equations similarly. Herein (D + k rx E) = o (V + k rx A) . 

7. CONCLUSIONS 

In the electromagnetic and gravitational fields, the field equations and their wave features can 
be rephrased with the algebra of quaternions and of octonions. The related conclusions include 
the wave equation, transverse wave, law of reflection, law of refraction, Fresnel formula, and total 
internal reflection, etc. 

With the algebra of quaternions, we can draw out the wave equation from the field equation 
in the gravitational field, and realize that the gravitational waves belong to transverse waves. 
Introducing the quaternion exponential function exp{i\a) , the wave vector and gravitational wave, 
g and bo , will produce two new wave components, g and b , respectively. Moreover, the g and 
b can be determined simultaneously, while the g and bo can be at the same time. 

In the electromagnetic field, making use of the algebra of octonions, we can deduce the wave 
equation, laws of reflection and refraction, Fresnel formula, and total internal reflection, and find 
that electromagnetic waves belong to the transverse wave in a vacuum. Bringing in the octonion 
exponential function exp(ila), the wave vector and electromagnetic wave components, Eo and Bo, 
will produce two new wave components, E and B , respectively. In contrast to the conventional 
electromagnetic theory with vector terminology, the research points out that the Eo and Bo can 
not be determined simultaneously via the wave vector. But the Eo and B can be measured at the 
same time, while the E and Bo can also be. 

Similar to the field strength, the force and torque etc may possess respectively the static quantity, 
alternating quantity, and independent wave etc. Besides the wave equation for field strength, there 
are also the wave equations for force and torque in the gravitational field, and that for force-like 
and torque-like in the electromagnetic field. In the octonion space, the mass continuity equation as 
well as the charge continuity equation are the gauge equations for Maxwell-like equations of torque 
wave and of torque-like wave respectively. The angular frequency oo of gravitational wave may be 
different to which of electromagnetic wave, force wave, force-like wave, torque wave, or torque-like 
wave etc. Similarly the speed vq of g/fo may be diverse from that of w e /vo or n e /vo etc, while the 
speed Vq of E/Vo may be distinct from that of W e /Vo or N e /Vo etc. 

It should be noted that the study for the deductions of wave equations in the electromagnetic 
fields and gravitational fields examined only some simple cases in quaternion spaces and octonion 
spaces. Despite its preliminary characteristics, this study can clearly indicate that the wave equation 
of electromagnetic field can be deduced from the field equation with the algebra of octonions, and 
obtain some inferences such as Fresnel formula, law of reflection, law of refraction, and total internal 
reflection etc in the vacuum far away from the field sources. For the future studies, the research 
will concentrate on only the electromagnetic waves transmitting in the conductors. 
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